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ON THE SPLASH SINGULARITY FOR THE FREE-SURFACE OF A 

NAVIER-STOKES FLUID 

DANIEL COUTAND AND STEVE SHKOLLER 


Abstract. In fluid dynamics, an interface splash singularity occurs when a locally smooth inter¬ 
face self-intersects in finite time. We prove that for d-dimensional flows, d = 2 or 3, the free-surface 
of a viscous water wave, modeled by the incompressible Navier-Stokes equations with moving free¬ 
boundary, has a finite-time splash singularity. In particular, we prove that given a sufficiently 
smooth initial boundary and divergence-free velocity field, the interface will self-intersect in finite 
time. 


1. Introduction 

1.1. The interface splash singularity. The fluid interface splash singularity was introduced by 
Castro, Cordoba, Feffernran, Gancedo, & Gomez-Serrano in [8j in the context of the one-phase water 
waves problem. As shown in Figure [TTT1 A splash singularity occurs when a fluid interface remains 
locally smooth but self-intersects in finite time. Using methods from complex analysis together 
with a clever transformation of the equations, Castro, Cordoba, Fefferman, Gancedo, & Gomez- 
Serrano [5j showed that a splash singularity occurs in finite time for the water waves equations. In 
Coutand & Slrkoller m, we showed the existence of a finite-time splash singularity for the one-phase 
incompressible Euler equations with free-boundary using a very different approach, founded upon 
an approximation of the self-intersecting fluid domain by a sequence of smooth fluid domains, each 
with non self-intersecting boundary. For one-phase flow, it is the vacuum state on one side of the 
interface which permits this finite-time interface self-intersection, and neither surface tension nor 
magnetic fields nor other inviscid regularizations of the interface change this fact mnn, and even 
stationary solutions, having a splash singularity, have been shown to exist (see Cordoba, Enciso, & 
Grubic 'TU]). 




Figure 1.1. The splash singularity at a point xq occurs when a locally smooth 
interface self-intersects in finite time t = T. 


On the other hand, for the two-phase incompressible Euler equations, wherein the moving interface 
is a vortex sheefi it was proven by Fefferman, Ionescu, & Lie [19] and Coutand & Shkoller m 
that a splash singularity cannot occur in finite-time while the interface remains locally smooth. In 
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particular, there is a fundamental difference in the behavior of the fluid interface when vacuum is 
replaced with fluid in the mathematical model. 

Since these results have been established for inviscid flows, it is natural to ask if splash singularities 
can occur for viscous flows modeled by the incompressible Navier-Stokes equations with a moving 
free-surface. Because the methods of constructing splash singularities for inviscid flows have relied 
on the ability to flow backward-in-time, a new strategy must be devised to study the parabolic 
Navier-Stokes equations. By using the change-of-variables employed in [8] together with stability 
estimates, Castro, Cordoba, Fefferman, Gancedo, & Gomez-Serrano in |Jj have shown the existence 
of finite-time splash singularities for the Navier-Stokes equations. Herein, we give a different proof 
which is amenable to any space dimension d > 2. 

1.2. The Eulerian description of the Navier-Stokes free-boundary problem. For 0 < t < T, 

the evolution of a d-dimensional (d = 2 or 3) one-phase, incompressible, viscous fluid with a moving 
free boundary is modeled by the incompressible Navier-Stokes equations: 


ut + u ■ V<r + Vp = vAu 

in 12(t), 

(1.1a) 

divu = 0 

in 12(f), 

(1.1b) 

v Def u • n — pn = 0 

on r(t), 

(1.1c) 

V(r(t)) = u- n 


(l.ld) 

U = Uq 

on 12(0), 

(l.le) 

12(0) = 12 0 . 


(l.lf) 


The open subset 12(f) C R d , d = 2 or 3, denotes the time-dependent volume occupied by the fluid, 
r(f) := 512(f) denotes the moving free-surface, V(F(f)) denotes normal velocity of r(f), and n(t) 
denotes the exterior unit normal vector to the free-surface T(t). The vector-field u = (ui,..,Ud) 

denotes the Eulerian velocity field, and p denotes the pressure function. We use the notation 

V = (di ,..., dd) to denote the gradient operator, and set Def u = Vu + Vri T , twice the symmetric 
part of the gradient of velocity. We have normalized the equations to have all physical constants 
equal to 1. 

The pressure p is a solution to the following Dirichlet problem: 

—Ap = u l ,ju\i in 12(f), (1.2a) 

p = n ■ \v Def u ■ n] on r(f) , (1.2b) 

so that given an initial domain 12 and an initial velocity field uq, the initial pressure is obtained as 
the solution of d) at t = 0. 

Definition 1.1. Given a locally smooth, time-dependent fluid interface or free-boundary, if there 
exists a time T < oo such that the interface T(T) self-intersects at a point while remaining locally 
smooth, we call this point of self-intersection at time T a “splash” singularity. 

We prove that there exist smooth initial data for the Navier-Stokes equations (0) for which such 
a splash singularity occurs in finite time. 

1.3. Statement of the Main Theorem. 

Theorem 1.1 (Finite-time splash singularity). There exist 

(1) open bounded C°°-class initial domains 12 C d = 2 or 3, with N denoting the unit normal 
vector field on <912, and 

(2) smooth divergence-free velocity fields uq satisfying the compatibility condition 

[Def uq ■ N] x N = 0 on <912, 

such that after a finite time T* > 0, the solution to the Navier-Stokes equations EJ has a splash 
singularity; that is, the interface T(T*) self-intersects. 
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In Theorem 110.11 we show that the geometry of such a splash singularity can be prescribed 
arbitrarily close (in the ff 3 norm) to any sufficiently smooth and prescribed self-intersecting domain. 


1.4. Prior results for the incompressible Navier-Stokes equations with moving free- 
surface. Local-in-time well-posedness of solutions to ED have been known since the pioneering 
work of Solonnikov |27H29j : his proof did not rely on energy estimates, but rather on Fourier- 
Laplace transform techniques, which required the use of exponentially weighted anisotropic Sobolev- 
Slobodeskii spaces with only fractional-order spatial derivatives for the analysis. Beale |5] proved 
local well-posedness in a similar functional framework, and Abels [T] established the existence theory 
in the L p Sobolev space framework. Well-posedness in energy spaces was established by Coutand 
& Shkoller in [12] for the case of surface tension on the free-boundary, and for Navier-Stokes fluid- 
structure interaction problems wherein a viscous fluid is coupled to an elastic solid, in PHH]. Guo 
& Tice [23] also used energy spaces for local well-posed for the case of zero surface tension. 

Beale [6] established global existence of solutions to ED for small perturbations of equilibrium. 
More recent small-data global existence and decay results (both with and without surface tension) 
can be found in ED, [26j i [25] i [20] , |4] . cind mm- Recent results on the limit of zero viscosity and 
the limit of zero surface tension can be found in ED, ED, and [32] . 

For the history of the well-posedness and singularity theory for the inviscid problem, we refer the 
reader to the introduction in [15] and HD 


1.5. Outline of the paper. In Section[2] we define our notation. In Section^ we define a sequence 
of domains f l e that we use as the initial data for the splash singularity, wherein the boundary T e 
of these domains is close to self-intersection with a distance e between two approaching portions of 
r e . We convert the Navier-Stokes equations to Lagrangian coordinates in Section [D thus fixing the 
domain. In Section [5] we present some preliminary lemmas which show that the constant appearing 
in elliptic estimates and the Sobolev embedding theorem is independent of e. In Section [G] we 
define the sequence of initial divergence-free velocity fields that are guaranteed to satisfy the single 
compatibility condition that we require, and whose norm is independent of e. Section[T]is devoted to 
the basic a priori estimates for the Navier-Stokes equations in Lagrangian coordinates; following our 
approach in [12], we establish estimates for velocity v G L 2 (0, T; H 3 (fl e )) fl C°([0, T]\H 2 (kl e )) which 
are independent of e. We then prove that the vertical component of velocity v(-t) at time t remains 
in an 0(t±) neighborhood of the vertical component of the initial velocity field. Using this fact, we 
prove the main theorem in Section (5[ we show that by choosing e appropriately, a finite-time splash 
singularity must occur at some time T* G (0, lOe). We consider a completely arbitrary geometry for 
a splash singularity in Section [9] by following our definition of a generalized splash domain from our 
previous work in [16j . This, then, allows us to show in Section m that we can construct a splash 
singularity for a geometry which is arbitrarily close in H 3 to any prescribed H 3 splash domain. 


2. Notation, local coordinates, and some preliminary results 


2.1. Notation for the gradient vector. Throughout the paper the symbol V will be used to 


denote the d-dimensional gradient vector V = ^ 


d d _d_ 

dxi ’ dx 2 ’ dx, 


d- 


2.2. Notation for partial differentiation and the Einstein summation convention. The 

kth partial derivative of F will be denoted by F,k = J^-. Repeated Latin indices i,j,k, etc., are 


summed from 1 to d, and repeated Greek indices a,/3,7, etc., 


example, F y , 


E d 
i =1 


and F\ a 1*00, p = £f =1 E„=r E 


are summed from 1 to d—1. 

2 aG ra/3 d(F_ 

0=1 dx a dxp * 


For 
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2.3. Tangential (or horizontal) derivatives. On each boundary chart Ui flO. for 1 < l < K , we 
let d denote the tangential derivative whose ath-component given by 

= ° er ‘= ( (V/oft) S 0(,r1 ' 

For functions defined directly on B + , d is simply the horizontal derivative d = ( d xi , ...,d Xdl ). 

2.4. Sobolev spaces. For integers k > 0 and a bounded domain U of K 3 , we define the Sobolev 
space H k (U) (if fc (/7;M 3 )) to be the completion of C°°(U) (C°°(U; ffi 3 )) in the norm 



for a multi-index a £ with the convention that |a| = 01+02 + 03 . When there is no possibility 
for confusion, we write || • ||fc for || • \\k,u- For real numbers s > 0, the Sobolev spaces H S (U) and 
the norms || • || Sj cr are defined by interpolation. We will write H S (U) instead of H S (U; K. d ) for 
vector-valued functions. 

2.5. Sobolev spaces on a surface T. For functions u £ H k (T), k > 0, we set 

ll u llfc,r = 

|a|<fc 

for a multi-index a £ For real s > 0, the Hilbert space H S (T) and the boundary norm | • | s is 
defined by interpolation. The negative-order Sobolev spaces H~ S (T) are defined via duality. That 
is, for real s > 0, H~ S (T ) = H S (T)'. 

2.6. The unit normal and tangent vectors. We let n(-, t) denote the outward unit normal vector 

to the moving boundary r(f). When t = 0, we let N e denote the outward unit normal to r e . For 
each a = 1 1 and x £ r e , T a (x) denotes an orthonormal basis of the (d— l)-dimensional 

tangent space to T e at the point x. 

3. The sequence of initial domains fT 

We shall use, as initial data, a sequence of domains, whose two-dimensional cross-section resembles 
a dinosaur neck arching over its body. 

3.1. The “dinosaur wave” domains. 

Definition 3.1 (The domain fl). Let H C R d , d = 2,3, be a smooth bounded domain (as shown on 
the left of Fiaure AS. 1\) with boundary T. We assume that there are three particular open subsets of 
Q as follows: 

(1) There exists an open subset w C H such that its boundary dui G T is a vertical circular 
cylinder of radius r and of length h > 0. 

(2) There exists an open subset C f1 which is the lower-half of an open ball of radius 1, 
located directly below the cylindrical region w, and in contact with the cylindrical region uJ. 
The “south pole” of lo+ is the point X + (see Fiaure [6J\) . 

(3) There exists an open subset w_ C LI directly below, at a distance 1, from the “south pole” 
X + of w +) such that the points with maximal vertical coordinate in dco- flT form a subset 
of the horizontal plane Xd = 0. 

(4) Coordinates are assigned to subsets of LI as follows: 

(a) The origin ofM. d is contained in C V fl {xd = 0}. 

(b) The point X +! the “south pole” of u> + , has the coordinates Xff = 0 for a = 1,..., d — 1 
and X d = 1 . 

(c) The top boundary of the hemisphere w+ is the set {( Xh,Xd ) £ R d : Xd = §, \xh\ < 1}- 

(d) The cylindrical region co is given by {( Xh,Xd ) £ R d : | < Xd < § + h, \xh\ < !}• 
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FIGURE 3.1. Left. The “dinosaur wave” domain Q with boundary T. Right. The sequence 
of “dinosaur waves” with boundary T e , e > 0, used as initial data for the Navier-Stokes splash 
singularity. In order to ensure that a splash occurs, the “dinosaur neck” ta e stretches downward 
so that there is a distance e between the two portions. The domains simply stretch the neck of 
the dinosaur, and are identical to fl away from the neck. 


Definition 3.2 (The initial domains IT). For 0 < eC 1 , let fl C R d , d = 2, 3, be a smooth bounded 
domain (as shown on the right of Figure 1,9.11) with boundary T e . We define the domain fl e to be the 
following modification of the domain fl: 

(1) There exists an open subset w e C fl e , which is a vertical dilation of the domain w, such that 
its boundary duj e fl T e is a vertical circular cylinder of radius r and of length h + 1 — e. 

(2) There exists an open subset u/j_ C IF which is the set translated vertically downward a 
distance 1 — e; hence, w!j_ is the lower-half of an open ball of radius 1, located directly below 
the cylindrical region w e , and in contact with the cylindrical region vj e . The “south pole” of 
o/j_ is the point Xf_. 

(3) There exists an open subset C fl e directly below, and a distance e, from the “south pole” 
XI of Lo e + , such that the points with maximal vertical coordinate in du>- flT form a subset 
of the horizontal plane Xd = 0. We assume that du>- fl F contains a d— 1-dimensional ball 
of radius sfe. 

(4) Coordinates are assigned to subsets of fl e as follows: 

(a) The origin of R d is contained in dco- C T fl {xd = 0}. 

(b) The point Xfi, the “south pole” of u/L, has the coordinates Xf = 0 for a = 1,..., d — 1 
and Xf = e. 

(c) The top boundary of the hemisphere is the set {( Xh,Xd ) el' 1 : Xd = e+ \xh\ < 
1}- 

(d) The cylindrical region uj e is given by {(xh,Xd) G R d : e+^ < Xd < e + ^+h, \xh\ < 1}- 

3.2. Local coordinate charts for fl and f Y. 

3.2.1. Local charts for O. We let s > 3 and 0<e<l. Let fl C denote a smooth open set, and 
let {U{\^L 1 denote an open covering of T = <9fl, such that for each l G {1,2,..., A'}, with 

B = 13(0,1), denoting the open ball of radius 1 centered at the origin and, 

B+ = BC{x d > 0}, 

B° = Bn{x d = 0}, 
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there exist C°° charts 81 which satisfy 

0;: B —> Ui is an C°° diffeomorphism, 

0i(b+) = Ui nn, 0 i(b°) = f/ ; nr, 


(3.1a) 

(3.1b) 


and det V 0/ = Ci for a constant Ci > 0. We assume these boundary charts can be split into three 
categories (each being non empty): 

• For 1 < l < K u 0i{B+) C uj. 

• For K\ + 1 < l < K 2 , 9 i(B + ) <£. uj and 9i(B + ) D u> + = 0. 

• For K 2 + 1 < l < K, di(B + ) <jt_ uj and di(B + ) D w+ ^ 0. 

We also assume that the images of any charts 6 i for K\ + 1 < l < K 2 does not intersect any of 
the images of the charts for K 2 + 1 < l < K. 

Next, for L > K. we let {Ui}\l k+1 denote a family of open sets contained in Q such that {Ui }[ =1 
is an open cover of LI and there exist smooth diffeomorpliisms 6 i : B —» Ui with det V 8 i equal to a 
constant Ci > 0 . 

Just as for the case of the boundary charts, we assume that these interior charts are split into 
three categories (each being non empty): 

• For K + 1 < l < Li, 8 i(B) C uj. 

• For L\ + 1 < l < L 2 , Qi(B) <jL uj and 9i(B + ) D w+ = 0. 

• For L 2 + 1 < l < L, 9i(B) <jL uj and 9i(B + ) fl 0. 


We assume that the union of the images of the charts 9i , for 1 < l < K\ and I\ +1 < l < Li contains 
the shortened cylindrical region uj= {( Xh,Xd ) : ^+e<Xd<^ + h — e, \xh\ < 1 }- 

We assume that the union of the images of the charts 6 i, for 1 < l < K\ and K + 1 < l < L\ 
contains the shortened cylindrical region 


w= {( x h ,x d ) G R d 


3 

2 


h h 

■7 : <x d < 


3 + h2 



h h 
3 + /i2’ 


\Xh I < 1} 


of length g 

We also assume that the union of the images of the charts 8 i, for K 2 +l < l < K and L 2 +l < l < L 
contains the complement in uj of the shortened cylindrical region 


UJ = {( Xh,Xd ) G : 


3 h ^ h 3 h ^ h 

2 + 2 3+T <Xd< 2 +h ~ 2 3+^’ 

\-\-h j 


\Xh\ < 1} 


of length , so that the complement is of length -§+ 7 ^. 

Finally, we assume the images of any of the charts 9i for L\ + 1 < l < L 2 does not intersect any 
of the images of the charts 81 for L 2 + 1 < l < L. 


3.2.2. Local charts for Q e . We next explain how this system of charts can be simply modified to 

describe fl £ using the following three steps: 

(1) For either 1 < l < K\ or K+l < l < L i, we define the vertically dilated chart (corresponding 
to a cylinder with length dilated from h to h + 1 — e) 

et=[e l ,e 2 , h + l h e (0 3 -^ + i + e ) . 

Note that 8 \ sends any point whose image by 9i was at the altitude § in uj (respectively 
| + h ) into a point of altitude \ + e (respectively | + h) in fl e . 

(2) For either K\ + 1 < l < K 2 or Li + 1 < l < L 2 , we set 8 ] = 8 i. 

(3) For either K 2 + 1 < l < K or L 2 + 1 < l < L, we set the translated in the vertical direction 
chart Of = 9i — (1 — e)e d - 

These charts describe f l e , and again detV0; is a strictly positive constant given by either Ci or 

h+l—e r' 

- ~ 
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3.2.3. Cut-off functions on charts covering D. Let {£;}/(-! denote a smooth partition of unity, sub¬ 
ordinate to the covering {Ui}^ =1 ; i.e., S 0 < & < 1 , and Y^iLi = 1 - 

We set Bj = B + for l = 1 and Bi = B for l = K + 1 For each l = 1 we set 

0 = o 9i, so that Q £ whenever the charts &i are smooth. 

3.2.4. Cut-off functions on charts covering f l e . We dehne the cut-off functions as follows: 

& 

Setting (i = £ £ o 0 ; e , we see that ||0||fc,Bi is bounded by a constant which is independent of e. 


4. The Lagrangian description of the Navier-Stokes free-boundary problem 


For e > 0, we let fT with boundary r e be given by Definition 13.21 and we transform the system 
ED into a system of equations set on this reference domain. To do so, we shall employ the 
Lagrangian coordinates. 

The Lagrangian flow map is the solution of the r) t (x,t) = u(r](x,t),t) for t > 0 with initial 

condition 77 ( 37 ,0) = 0. Since divu = 0, it follows that detV ?7 = 1. For each instant of time t for 
which the flow is well-defined, we have 

r](-,t) : D £ —> Q(t) is a diffeomorphism; 

furthermore, thanks to (I1.1H L 

r(t) = ?? (r e ,t). 

Notationally, we keep the dependence on e > 0 implicit, except for the initial domain and boundary. 
Next, we define 


v = uo rj (Lagrangian velocity), 
q = p o 77 (Lagrangian pressure), 

A = [Vr ;]^ 1 (inverse of the deformation tensor), 
g a) 3 = r], a -T],p a, f3 = 1 ,.., d — 1 (induced metric on T), 

0 = det(< 7 a| a). 

We also define the Lagrangian analogue of some of the fundamental differential operators present in 
this equation: 

div,, v — (div u) o 77 = v l ,j Aj , 
curl,, v = (curlit) o 77 or [curl,, v]i = £ijkV k , r A £ , 

Def,, v = (Def u ) o 77 or [Def,, v]j = v l , r A J + v\ r A\ , 

A v v = (A u) o rj= (. A 3 r A k v, k ),j ■ 


The Lagrangian version of equations ED is given on the fixed reference domain f l e by 

r t 


r](-,t) = e+ / v(-,s)ds infFx[0,T], (4.1a) 

Jo 

vt + A T Vq = uA v v in x (0, T ], (4-lb) 

div,, v = 0 in fT x [0, T ], (4.1c) 

uDef,, v ■ n — qn = 0 on F e x [0, L], (4. Id) 

( 77 , v) = (e, u 0 ) in fT x {f = 0} , (4.1e) 


where e(x) = x denotes the identity map on Q, and where we write n for 77 .( 77 ) in the Lagrangian 
description; in particular, the unit normal vector n at the point rj(x, t ) can be expressed in terms of 
the cofactor matrix A and the time t = 0 normal vector N e as 


n = A T N e /\A T N £ \. 
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Due to dnt), 

A n v = div, ; Def^ v , 

so that (|4.1H 1 can be viewed as the natural boundary condition. The variables 77 , v, and q have an 
a priori dependence on e > 0 , but we do not explicitly write this. 

Local-in-time existence and uniqueness of solutions to ed have been known since the pioneering 
work of Solonnikov [27] . We shall establish a priori estimates for ED with the initial domain 12 e 
and with divergence-free initial velocity fields satisfying the single compatibility condition 

[Def u% ■ N e ] ■ = 0 on F e , (4.2) 

where N e denotes the outward unit normal to T e and a = 1,d — 1, denotes the d— 1 tangent 
vectors to T e . 

We will show that both the a priori estimates and the time of existence for solutions are indepen¬ 
dent of the distance e > 0 between the falling dinosaur head X^_ and the flat trough du>- D {xa = 0 } 
(see Figure [Tip . To do so, we shall rely on some basic lemmas that provide us constants which are 
independent of e. 


5. Elliptic and Sobolev constants are independent of e 
We consider the following linear Stokes problem 


A u + Vp = / 

in 12 e , 

(5.1a) 

div u = <p 

in 12 e , 

(5.1b) 

u = g 

on r e , 

(5.1c) 


Lemma 5.1 (Estimates for the Stokes problem on IT). Suppose that for integers k >2 , f £ 
H k ~ 2 (TL e ), (f> £ iL fe_1 (12 e ), and g £ H k ~ 1 ^ 2 ( T e ), and f Qe cf>(x)dx = J re g-NdS. Then, there ex¬ 
ists a unique solution u £ H k (tt e ) and p £ iJ fe_1 (12 e )/M to the Stokes problem \5.1\) . Moreover, 
there is a constant C depending only on 12, but independent of e > 0, such that 

IMU + IblU-i < c (II/IU -2 + II 0 IIfc —1 + lflU- 1 / 2 ) • (5-2) 

Proof. The estimate (15.211 is well-known on the domain D; see, for example, [2]. This estimate on 
the sequence of domains 12 e follows by localization using the charts d\ given in Section 13.21 Since 
the charts 9\, are modified from the charts 9i by a vertical dilation with lower and upper bound that 
is uniform in e, the constant for the elliptic estimate in each chart is independent of e > 0. □ 

Lemma 5.2 (Sobolev constant on 12 e ). Independent of e, there exists a constant C > 0 which depends 
only on the domain 12, such that 

max |u(a;)| < C||u|| s ,n e VuSiL s (12 e ), s > d/ 2. 

Proof. The constant is determined by the radius r of the smallest ball B(x, r) for x £ 12 e , such that 
B(x, 7') C 12 e . By Definition 13.21 of the domains 12 e , r does not depend on e, and hence the Sobolev 
constant C only depends on 12. □ 

Lemma 5.3 (Trace theorem on 12 e ). Independent of e, there exists a constant C > 0 which depends 
only on the domain 12, such that for s £ (|, 3] 

NL-i,r« < C1M| a , n . Vu£H s ( 12 e ). 

Proof. From the standard trace theorem in B + , we have the existence of a constant C > 0 such that 
for any boundary chart, 


H«o0fL_i iBo < C\\uo9l\\ StB+ Vu£H°( 12 e ). 
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Now, since 9 is either a chart 9i for the domain U or a vertical dilation of such a chart with a 
uniform bounded from below and above as is made precise in Section 13.21 this implies that by the 
chain rule, 


IMIs—±,0f(.B o ) - C\\u\\ s fie( B +) 


VueH s (n e ). 


Since T e is the union of all 91 (Bo), 1 < l < K, the above inequality implies the result. 


□ 


6. The sequence of initial velocity fields Uq 

6.1. Constructing the sequence of initial velocity fields Uq. As described in Definition 13.21 
near the intended splash (or self-intersection) point, the open set f l e consists of two sets: the upper 
set and the lower set w_ whose boundary contains the flat “dinosaur belly” at Xd = 0, as shown 
in Figure IfTTl We We let denote the point which has the smallest vertical coordinate in du>+. 
Directly below, we let X- be the point in did- fl {xd = 0} with the same horizontal coordinate as 
XI. Without loss of generality, we set X_ to be the origin of R d . 



FIGURE 6.1. In a neighborhood of the intended splash point, we suppose that iV consists of 
two sets: the upper set cjf, and the lower set uj— containing the horizontally flat “dinosaur belly.” 
The point X+ is at a distance e from the set uj— and the point X— is assumed to be the origin in 
K d . 


We choose a smooth function 6g G C°°(r e ) such that &§ = — 1 in a small neighborhood of XI on 
du> e + , = 0 on duj-, = 0 on <9w e fl r e , f re bg dS = 0, and satisfying the estimate 

||froll2.5,r* < m 0 < oo , (6.1) 

where mo does not depend on e. 

We define the initial velocity field Uq at t = 0 as the solution to the following Stokes problem: 


—Auo + V?o = 0 in , (6.2a) 

diviig = 0 in Q e , (6.2b) 

[Def Uq ■ A e ] • = 0 on T e , (6.2c) 

u e 0 ■ N e = b e on r e , (6.2d) 


with N e denoting the outward unit normal to T e and a = 1,2 denoting an orthonormal basis of 
the tangent space to T e (if the dimension d = 2, then there is only one tangent vector). Using the 
regularity theory of this elliptic system (see, for example, [3U] or [3] and references therein), together 
with the proof of Lemma 15.11 for a constant independent of e > 0, 


||woll3,fi e < C'||6 e ||2.5,r« < Cmo. 


(6.3) 


The boundary condition (16.2t b ensures that Uq satisfies (14.21) . 
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6.2. The initial pressure function p q. The initial pressure function pg at t = 0 then satisfies 

-ApS = {u e 0 )\j (uoYn in ST, (6.4a) 

Pq = Nq ■ [v Def Uq • Nq] on T 6 , (6.4b) 

so that using the same proof as that of Lemma 15.11 we have the following e-independent elliptic 
estimate: 

IIP5ll2.fl- < C [Klkn. + Ikells.ne] < C7P(mo), (6.5) 

where we use V denote denote a generic polynomial function that depends only on S l. 


7. A PRIORI ESTIMATES 

Let ST denote the dinosaur domain shown in Figure 13.11 and let 9i denote the system of local 
charts for ST as defined in m- By denoting rji = p o 8i we see that 


rji(t) : B + —► Q(t) for l = 

We set Vi = uory, qi = pory and Ai = [Dp/] , fi = Ci (where Ci > 0 is a constant, and ai = JiAi. 


The unit normal ni is defined as g 2 x plL if d = 3 and by g 2 if d = 2. 


It follows that for l = 1, K, 

Vi ( t ) =0i+ f vi 

Jo 

in B + x [0, T ], 

(7.1a) 

dm +AjX 7 qi= A m vi 

in B+ x (0, T \, 

(7.1b) 

O 

II 

S? 

> 

‘^3 

in B+ x [0, T], 

(7.1c) 

v Def^, Vi ■ ni - qini = 0 

on B° x [0, T], 

(7-ld) 

(VhVi) = ( 0 h uo 0 81) 

in B + x {t = 0} , 

(7-le) 


where we have set v = 1. 


Definition 7.1 (Higher-order energy function). For each t £ [0,T], we define the higher-order 
energy function 


E e {t) = 1 + |W-,i)||3,n« + IK:*)ll!,fi‘ + f IK 

Jo 


s )ll3,O e ^ s 


IKKH^ds 


lkM)llo,n« 




We then set Mo = V(E e ( 0)) where V denotes a generic polynomial whose coefficients depend only 
on SI. The constant Mo is then equal to V(mo), a polynomial function of the constant mo introduced 
in \6.3\) . 


Theorem 7.1. Assuming that T(t) does not self-intersect, independent of e > 0, there exists a time 
T > 0 and a constant C > 0 such that the solution 

v £ C([0,T],ff 2 (ST))nL 2 (0,T;f7 3 (ST)), q £ L 2 (0, T; tf 2 (ST)) 

to HIP satisfies the a priori estimate: 

max E e (t) < C Mq . (7.2) 

te[o,T] 

Proof. The proof will proceed in five steps. 


Step 1. Estimates for Vp and A. Using (|7. lh l . we see that 


||V?/(-,f) -Id|| 2 ,o. < 


Vv(-, s)ds 


< Vt sup \JE e (t ) . 
2 ,o e «e[o,t] 


(7.3) 
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Thanks to Lemma 15.21 there exists a constant C > 0, independent of e, such that 

l|Vry(-,t) -Id|| L =o (ne) < cVt sup ^E*(t). 

se[o,t] 

Since det V 77 = 1, the matrix A is simply the cofactor matrix of V 77 : 


(7.4) 


A = 


~V,2 

v,i 


for d = 2, and A = 


Vi2 X?7,3 

??,3 >< 77,1 ford = 3, (7.5) 

Jla x?y,2_ 

where each row is a vector, and for a 2-vector x = ( Xi,X 2 ), x ± = (— X 2 ,Xi). 

We make the following basic assumption, that we shall verify below in Step 5: for a constant 
0 < d -C 1, we suppose that t £ [0, T] and that T is chosen sufficiently small so that 

(7.6) 


sup ||Vry(-,t) — Id ||L=o ( f 2 e) < d 10 . 
te[o,T] 

It follows from (EH), that since ||A(-, t) - Id ||z,~ ( n« } < fo Pt(u s)\\ L ^ {n ^ds, 

sup ||A(-,t)-Id lUoop,.) + \\AA T (-,t) — Id ||z,°°(n e ) < $ ■ 
te[o,T] 


(7.7) 


Step 2. Boundary regularity. We begin by considering a single boundary chart 9i : B + —> 12 (f). 
Let Q denote the smooth cut-off function defined in Section l5.2.4l Using equation (17.1b ). we compute 
the following L 2 {B + ) inner-product: 

(0 d 2 [d t vi - A v v + Aj V 9; ] , Cid\) L2{B+) =0. (7.8) 

To simplify the notation, we fix l £ {1,..., K} and drop the subscript. The chart 9i was defined so 
that det Vd; = C; for a constant Ci > 0. Then m can be written as be written as 


'B+ 


C, 2 d 2 v\ d 2 v l dx — / ( 2 d 2 [A^A J s v\j],k d 2 v l dx + / £ 2 d 2 [A k q\, k 8 2 v l dx = 0. (7.9) 


lB+ 


0=~Kd 2 v(ml B+ + B 2 [A k s Aiv i , j ]8 2 [C 2 v%dx+ I 8 2 [A k q]8 2 [C, 2 v% k dx (7.10) 


Integration-by-parts with respect to Xk shows that 
1 d 

--IKouuUJIio, 

^ J B+ J B+ 

where we have used the boundary condition (17.m i to show that the boundary integral vanishes. 
Using <P' fc to denote the Kronecker delta function, we write (17.101) as 
1 d 


2 dt 


IIW,*)llo, B+ + IIC^Vi;(t)||^ + =- / d 2 [A k q]d% 2 v\ k dx 


lB+ 


'B + 


8 2 [(A k s Ai - S k i)v\i}8 2 [C 2 v%kdx- [ [d 2 v\ k (8 2 C 2 v i + 28( 2 8v l ),k +U8 V] dx. 

Jb+ 

(7.11) 


We integrate (17.111) over the time interval [0,T]: 


-||Ca 2 i;(-,t)||g iB+ + / W(d 2 vmlB + <M 0 +i 1+ i 2 +i 3 


where 


(7.12) 


Ti = 


I 2 = 


I 3 = 


/ 0 Jb+ 

f 

/ 0 J B+ 

n 

/ 0 J B+ 


8 [A k q] d 2 [C v l ] lk | dxdt , 

1 8 2 [(A k s Ai - ] 8 2 [C 2 v% k I dxdt , 


d 2 v l , k [d 2 ( 2 v' t + 2i d( 2 dv l ],k +£, k d 2 v l \ dxdt. 
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Using the Sobolev embedding theorem and Lemma 15.21 We estimate I\ 


h < 


/ 0 JB+ 


\d q\ \A k d 2 v\ k | dxdt + / ||g|| 2 , £ |m| 2 ,n'|Mkn*di 


A 


+ / |M|l.5,e|m|2,n«|M|3 > fi e di • 


To estimate the integral If, we use (EHO to write 

n ,i Ak _ Ak n ,i \k n ,i \k n A 

V ikot.fi -A-i — -^-i iotfi V ifi V ikot iot V ikfi , 

so that the term with three derivatives on v is converted to a term with three derivatives on 77 plus 
lower-order terms. It follows that for <5 > 0, and a constant C$ (which blows-up as <5 — > 0), 

li < 5 [ \\q\\ln.dt + C S TP( sup E\t )). 

J o te[o,T] 

The integral if is estimated in the same way. For the integral If we use linear interpolation to 
estimate the norm f Q ||q|| i. 5 , e : 

Zf < S [ T |Ml!, n <dt + S [ T \\q\\l n .dt + C S TP( sup E\i)). 

Jo Jo te[o,T] 

It follows that 

Ii < M 0 + C S TP{ sup E e (t)) + S sup E e (t). (7.13) 

te[o,T] te[o,r] 

Next, for the integral I 2 , 

r T r rT 

I 2 < / 


to J B + 


I (A*A* - 6 kj )d 2 v i , j d 2 [C, 2 v% k | dxdt + 2 f f \d(A k A> s - 5 k: >)dv\ 0 d 2 [<?v% k \ dxdt 

Jo Jb+ 


A 


+ 


to Jb+ 


I B 2 (A k Al - 5 k i)v\ 0 d 2 [( 2 v% I dxdt . 


Using m and choosing 1 9 < 6, 

I 2 a < CsTP( sup E e (t)) + 5 sup E e (t). 
te[o,T] te[o,T] 

In the same way as above, we again use Lemma 15.21 together with linear interpolation for term If, 
to see that 

I 2 <M 0 + C S TP{ sup E e (t)) + 5 sup E e (t). (7.14) 

te[o,T] te[o,T] 

The integral I 3 is straightforward and also satisfies 


I3 < Mo + CsTP( sup E e (t)) + CS sup E e (t). 

te[o,T] te[o,T] 


(7.15) 


Summing over all of the boundary charts l = 1 in (17.121) . the inequalities (17.131) (17.151) 

together with the trace theorem, Lemma 15.31 show that 

[ IN", t)\\l. 5 ,r* < Mo + CgTP( sup E e (t)) + 5 sup E e (t) 

Jo ’ tefo.Tl tefo.Tl 


(7.16) 
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Step 3. Estimates for the time-differentiated problem. We consider the time-differentiated 
version of m which we write as the following system: 


V t =V 


in n e x [0, T ], 

(7.17a) 

+ A T Wq t = -AjS7q + [d t (AiA k 8 )v, k ] tj 

in IT x (0, T} , 

(7.17b) 

div^ v t = -v\j d t Aj 


in n e x [0, T ], 

(7.17c) 

> ■ n — qn] = 0 


on r e x [0, T ], 

(7.17d) 

cT 

g 

£ 

II 

£ 

O 

e 

i—* n\ 


in x {t = 0} , 

(7.17e) 


where u\ = Aitg — Vpjy with u e a defined in (16.21) and defined in (16.41) : therefore, independently of 
e > 0, 

||«il|o,n. <V{m 0 ). (7.18) 

We define the space of div^-free vectors fields on f l £ as 

V(t) = {cj,GH 1 (^-R d ) : div„ ( ., t) <j> = 0} . 

Taking the L 2 (fl e ) inner-product of equation (17.17b ) with a test function </> € V(t), we have that 


/ v tt ■ <t>dx + / d t [A k s Alv\j\(j)\ k dx = / qd t A!?4> l , k dx V0 G V(t). 

Next, we define a vector field w satisfying 

div,, w = — v l ,j dtA\ in f, 
w = <f>{t)n on T e , 


(7.19) 


(7.20a) 

(7.20b) 


where = — f Qe —v\j dtAldx/\T e \. A solution w can be found by solving a Stokes-type problem, 
and according to the proof of Lemma 3.2 in HU. for integers k > 1, 

IM-,t)IU,n* <c(^\v\ 3 (-A)d t Al(-,t)\\ k -i&. + \\(t>{t)n\\ k _ l/2X ?j , (7.21) 

where the constant C is independent of e by Lemma PTTTl It follows from (17.211) and (17.51) that 

sup \\w{-,t)\\ 2 ,n'+ [ |K-,t)||! in . <TP( sup E e (t)). (7.22) 

tg[0,T] JO tgfO.Tl 


Similarly, 


div^ w t = - (w\j d t Aj + d t (v\j d t Aj )) in fl e , 
w t = {(j>tn) t 


on r e . 


(7.23a) 

(7.23b) 


and 


so that 


IK||i,ne < C (|| w l ,j d t Aj + dt(v l ,j dtAj) ||o,n e + ||(</ > t”-)t||i/ 2 ,r' 


/ IKIIi,n« < TP( sup E e (t)). 

Jo tefo.Tl 


(7.24) 


te[o,T] 

Now, because of (17.20b ,). v t — w£V(t), and we are allowed to set (j) = v t — w in (17.191) . We find that 
+ f dt[AgA J s v l ,j ] v\,k dx = f v tt -wdx+ f d t {A k a A : > s v t , j )w t , k dx 

+ f qd t Ai [v l t , k +w\ k ] dx. 

Jq 










14 


D. COUTAND AND S. SHKOLLER 


and hence for t £ (0, T), 


Ji 


IM->*)llo,o 6 + / HV« t ||§, n .da = ^IMIo,n e ~ / / - S kj ]v % t ,j v\, k dxds 

Jo z Jo J n e 


r* 

/0 ./fi' 


5t[^^]u 8 ,j^,feda:ds- 


• wdxds + 


o 


o Jiv 


dt[A k A J a v l ,j ] w’,fc dxds 


Ji 
1 k f 


J 3 


Ji 


o Jn 


qd t A k [v\, k +w\ k ] dxds . 


Js 


For S > 0 and using (17.711 with d < S, we see that 


|Ji| < S sup E e (t). 
te[o,T] 


(7.25) 


Next, according to m the components of A are either linear ( d = 2) or quadratic (d = 3) with 
respect to the components of V?;; hence, dtA behaves like Vu for d = 2 and like Vr/ Vu for d = 3. 
We consider the more difficult case that d = 3 in which case dt(AA T ) behaves like V?/ V?/ V// Vu. It 
follows by the Cauchy-Young inequality that for S > 0, we have that 


\Ji\ < Mo + TV( sup E e (t))+S sup E e (t). 

te[o,T] *e[o,T] 


(7.26) 


To estimate Js, we integrate-by-parts in time: 


\Jz\ < / / \vfW t \dxds 

Jo J n e 


i>t ■ wdx 


< 


[ f \v t ■ w t \dxds + M 0 + f \v t {-,t)w(-, 0) |<ia; + f f w t (-,s)ds 

Jo Jit* Jn e Jn e Jo 


dx 


<M 0 + TV( sup E e (t))+6 sup E e (t ), 

te[o,T] te[o,r] 


(7.27) 


the last inequality following from the Cauchy-Young inequality and the estimates (17.211) and (17.241) . 
The integrals Ji and Jo (using (17.221) and (17.241) ') are estimated in the same way as Ji so that 


\Ja\ + \J 5 \<M 0 + TV( sup E e (t)) + 5 sup E e (t). 

te[o,T] te[o,T] 


(7.28) 


Combining the estimates (17.251) (17.2811 . we find that 


sup ||t>t(-,i)||o,fie + [ \\vt\\l,nedt < M 0 +TV{ sup E e (t)) + C6 sup E e (t) . (7.29) 

t€[ 0 ,T] Jo telO.Tl tefo.ri 


te[o,T] 


te[o,T] 


Step 4. Regularity for the velocity and pressure. Next, we write equation (14.1b 'l as 

—Av + Vij = div[(dd T — Id)Vu] — ( A T — Id)Vg — v t in f Y x (0, T] , (7.30a) 

divu = — (A? — in Q e x [0,T], (7.30b) 

v € L 2 (0,T;R 2 - 5 (T e ) (7.30c) 

The two inequalities (17.161) and (17.291) . together with the Stokes regularity given in Lemma l5Jl 
show that v € i°°([0, T\; H 2 (fl e )) D L 2 (0, T; H 3 (tt € )) and satisfies 


sup \\v(;t)\\ln* 

te[o,T] 


+ / IMl3,n eC ^+ 


| l n .dt< M 0 +TV( sup E e (t))+CS sup E e {t) . (7.31) 


te[o,T] 


te[o,T] 
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By choosing S > 0 sufficiently small, we obtain that 

sup E e {t) < M 0 + TV{ sup E e (t)), (7.32) 

te[o,T] te[o,T] 

for a constant M$ and a polynomial function V which are both independent of e. 

From the estimate (17.311) . v G L 2 (0, T; 77 3 (fF)), and the estimate (17.291) . v t G L 2 (0,T; i7 1 (fl e )). 
Using the partition of unity functions Q defined in Step 2 above, we then see that for each chart 
Qv G L 2 (0, T; H 3 (Bi)) where Bi = for l = and Bi = B for l = K + 1 Similarly, 

Qv t G L 2 (0, T; H 1 (Bi)). It is then standard that Qv G C°([0, T]; H 2 (Bi)) : and hence by summing 
over l = 1, v G C°([0, T]; 7J 2 (fl £ )). 

Since the pressure satisfies the elliptic system: 

= u*, r A^v 3 ',* A* in x (0, T], 
q = n- [Def,,!; • n) on T e x [0, T ], 

we then infer that q G C^QO, T]; iJ 1 (U e )). Then, using the momentum equation (17.30al) . it follows 
that Vt G C'°([0, T]; L 2 (f2 e )). 

This then shows that E e (t) is a continuous function of time. Following Section 9 in [14], from 
(17.321) . we now may choose T > 0 sufficiently small and independent of e, such that 

sup E e (t)<2M 0 . (7.33) 

te[o,T] 

Step 5. Verifying the basic assumption (17.61) . Having established (17.331) on [0,T] with T 
independent of e, for any e > 0, we may now use the formula (17731) to choose T even smaller if 
necessary to ensure that (17.61) holds. This concludes the proof. □ 

We now establish a more quantitative estimate in order to assess the continuity of d 2 v(t,-) in 
i 2 (U e ). 

Proposition 7.1. For all t G [0,T], 

mtut\\d 2 (v e (-,t)-u e 0 )\\l nc + f \\d 2 (v t (-,s)-u e 0 )\\ 2 in eds<t 1/2 V(M 0 ). (7.34) 

sen* j 0 

Proof. We write v(t) = v(-,t) and again set viscosity v = 1. The difference v(t) — Uq satisfies the 
equation 

(v - u e 0 ) t - A v (v - Uq) + A T Vq = A v u e 0 . 

Following Step 2 in the proof of Theorem 17.11 and once again localize to a boundary chart 6 i , 
l = 1 with det = Ci and with cut-off functions (i, we obtain that 

° = \jt KB2[v{t) ~ + J B+ ^ A ‘ A ^ V - 1 • ^[C 2 (- - u o)U dx 

+ [ B 2 [A k l q\5 2 {f 2 v l ) lk dx+ [ d 2 [A k s Alu e 0 , j ]-d 2 [C 2 (v-u € 0 )], k dx, (7.35) 

JB+ JB+ 

where we have dropped the explicit chart dependence on l and where again, the boundary integral 
terms have vanished due to (14.11 11. We integrate (17.351) over the time interval [0,T]: 

||C<5 2 K*) - Uo]IIo,b+ + / IIC<5 2 b(f) - Uo]lll,B+ < |/Ci| + |/C 2 | + \)C 3 \ + |JC 4 |, 

Jo 
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where we are writing Uq for Uq o Oi , and where 

pT p 

ic 1 = / 


/ 0 JB+ 
,-T r 


AC:2 = 


Jo J B + 

ic,= r [ « 2 


/ 0 ./B+ 

pT 


B 2 [A k q] 3 2 {C, 2 {v - u e o y], k dxdt , 

-a 2 [c 2 (^-«o)]>fc rfa;di . 

a 2 (w - «o)\ fe [ [<9 2 C 2 0 - '«o) 1 + 2 d( 2 d(v - u e o y], k +C 2 ,fe B 2 v l ]dxdt , 


K a = 


d 2 U k s A{ul 0 ]-d 2 [( 2 (v~ul)Udxdt. 


) 0 J B + 


We write 


ACi < 


/o Js+ 


<9 [A*g] <9 2 [£ V], fc cfedt + 


/o Js+ 


a 2 [A^]9 2 [c 2 <], fe 


dxdt . 


/c? /c; 

By (16.31) and (17.21) . we see that 

|ACj| < Vtp(m 0 ). 

For the integral AC“, we focus on the integrand that arises when B 2 acts on both q and v l , k , for all other 
derivative combinations immediately give an integral bound of VtV(Mq). Using the Lagrangian 
divergence-free condition (14.1b ). 


r T f 


f T f 

/ / C 2 B 2 qA k 3 2 v\ k dxdt 

< 

/ / C 2 B 2 qB 2 A k v\ k dxdt 

Jo J B+ 


Jo J B+ 


/ / C, d q 3A^dv l , k dxdt 

'0 J B+ 

An application of the Cauchy-Young inequality together with the Sobolev embedding theorem, shows 
that 

|ACJ| < Vtv{m 0 ) . 

For the integral 1C 2 , we consider the case that B 2 acts on (A k A J s — S k l), all other terms immediately 
giving the desired bound. Using (17.31) and (17.51) . ||AA T —Id || L oo( B +) < \/TV(M), so that with m , 

\ic 2 \ < Vrr(Mo). 

The integral 1C 3 and 1C 4 are easily estimated using the Cauchy-Young inequality, the Sobolev em¬ 
bedding theorem, and (17721) . We have thus established that 

Iica 2 [u,(f) - 4 ]|| 2 , s+ + [ T \\cB 2 [ Vl (t) - U * 0 o e l ]\\l B+ < Vtv(Mo) . 

Jo 

Summing over l = 1 then concludes the proof. □ 

8 . Proof of the Main Theorem 

Using the Lagrangian divergence condition (14. Ib ). we have that divw = —(Aj—dj)v l ,j, which we 
write as divu = —(A — Id) : Vv. Then, since divug =j for all t € [0,T], 

||ddiv(u — tto)llo,n 6 < \\B(A - Id) Vu||^ ne + ||(A - Id) <VfV{M 0 ). ( 8 . 1 ) 

Using EB together with (17.341) . the normal trace theorem (see, for example, (A. 6 ) in [16]) shows 
that B 2 (v — Uq) ■ N e £ C([ 0 , T; (r e )) and 

\\d 2 (v-u^-N e \\ 2 _ 1 / 2 ^<VTV(M 0 ), 
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SO that 

||(«-«S)-JV e ||?. 5ir . < Vtv(Mo), 

and hence by Lemma 15.21 

max.\\(v(x,t)— Uq) ■ N e \\<T*V(Mo) Vi£[0,T]. (8.2) 

a?Gr e 

Next, we consider the motion of the points X f \ and given in Section ROI (see Figure 1(01) . 
Recall that the unit normal N e at both the points XI = (0,0, e) and X _ = (0,0,0) is vertical, so 
by definition of Uq, we have that 

Uq(X^_) • N e = — 1 Uq(X-) ■ N e = 0 , and \X+ — X- \ = e . 

Using Theorem l7.ll we choose e so small that lOe < T, where [0,T] is the time interval of existence 
which is independent of e, and we consider the vertical displacement of the falling particle X+. Since 
X • e d = e, and 

ri(X^,t) ■ e d = e+ [ v d (X^_,s)ds, 

Jo 

for t = lOe, we have from (18.21) that 

r) d (X+, lOe) < —8e . 

Next, let Z denote any point on duj- D {x d = 0}. Since Uq(Z) ■ N e = 0 and rj(Z, lOe) = v(Z, s)ds, 
according to <E2D, 

r/(Z, lOe) • e d > —cei , c = 10* V(M 0 ). 

We then choose e > 0 sufficiently small so that cei < 8e. It follows that 

v(X+, 10e) ■ e d < ri(Z, lOe) • e d ■ (8.3) 

We next consider the horizontal displacement of the particle X f + and any particle Z on duj-C\{xd = 
0} x [0,10e]. From the estimate (17.331) . for all time t £ [0, lOe], ||u(-, f)||,L°o(Q) < V(M 0 ). 

Therefore, for any t £ [0, lOe] and for a = 1,..., d— 1, 

\r) a {X%,t)\ < 10eiP(Mo) and \r] a (Z,t) - Z a \ < 10 eV(M 0 ), 

showing that the distance between the projection of the surface r](d w_ D {x d = 0}, t) onto the plane 
x d = 0 and the set <9w_ D {x d = 0} is 0(e). Since by Definition 13.21 the set doj- D {x d = 0} contains 
a d— 1-dimensional ball of radius ^fe centered at the origin, we see that by choosing e sufficiently 
small the vertical line passing through rj(X+,t) must intersect the surface r](du}- D {x d = 0},f) for 
any t £ [0, lOe]. Now, since at t = 0, X+ is directly (vertically) above cU>_ D {x d = 0}, and at 
t = lOe, from (18.31) . r](X+, lOe) is (vertically) below r/(du}- D {x d = 0}, lOe), then by continuity there 
necessarily exists a time 0 < T* < lOe at which T*) = r](Z , T*) for some Z £ D {x d = 0}. 

This concludes the proof of the main theorem. 

9. The case of a general self-intersection splash geometry 

We now show how the analysis presented in the previous sections for the case of the “dinosaur 
wave” initial domain can be used to establish the existence of a splash singularity in a finite time T* 
for any domain whose boundary is arbitrarily close (in the i7 3 -norm) to any given self-intersecting 
surface of class H 3 . This generalization requires the geometric constructions that we introduced in 
our previous work m, coupled with a very minor adaptation of the analysis of the previous sections. 
We begin with the definition of the splash domain that we gave in m■ 
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9.1. The definition of the splash domain. 

(1) We suppose that xq £ T := is the unique boundary self-intersection point, i.e., f l s is 
locally on each side of the tangent plane to dil s = T s at xq. For all other boundary points, 
the domain is locally on one side of its boundary. Without loss of generality, we suppose 
that the tangent plane at Xq is the horizontal plane X3 — (£ 0)3 = 0. 

(2) We let Uq denote an open neighborhood of Xq in R 3 4 5 , and then choose an additional L open 
sets {Ui}f =1 such that the collection {f//}^L 0 is an open cover of T s , and {Ui}^L 0 is an open 
cover of and such that there exists a sufficiently small open subset wCl/o containing xq 
with the property that 

aJD Ui = ® for all l = l,...,L. 

We set 

Uq = U 0 n n s n (a; 3 > ( x 0 ) 3 } and Uq = U 0 n f2 s n (a; 3 < ( 3 : 0 ) 3 } ■ 

Additionally, we assume that Uq D fl s D {a ; 3 = ( 2 : 0 ) 3 } = {x’o}, which implies in particular 
that Uq and Uq are connected. See Figure 9.1. 



Figure 9.1. Splash domain Q s , and the collection of open set {Uo, C/ 1 , C/ 2 , Uk} 

covering T. 

(3) For each l £ {1,..., K }, there exists an H 3 -class diffeomorphism 9i satisfying 

0i : B := 5(0,1) -A U t 
u l nn s = d l (B+) and ^nr, = 0,(B o ), 

where 

B + = {(aq ,x 2 ,x 3 ) £ B : x 3 > 0} , 

B° = {{x 1 ,x 2 ,x 3 ) £ B : x 3 = 0} . 

(4) For L > K, let { U,}^ K+l denote a family of open sets contained in Q s such that (L}}^_ 0 is 
an open cover of f2 s , and for l £ {A' + 1,..., A}, 0; : B —> Ui is an H 3 diffeormorphism. 

(5) To the open set Uq we associate two A 3 -class diffeomorphisms 9 + and 6- of B onto Uq with 
the following properties: 
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e + (B+) = u+, e_(B+) = u^, 

9+(B°) = u+n r s , e_(B°) = u^n r s , 

such that 

{xo } = 9 + (B°)ne.(B°), 

and 

0 + (O) = MO) = *o. 

We further assume that 

6 ±(B+ n B{ 0,1/2)) n 9i{B+) = 0 for 1 = 1 

and 

9 ± (B+nB(0 1 l/2))<l9i(B) = 0 for l = K +1, ..., L . 

Definition 9.1 (Splash domain 12 s ). We say that Cl s is a splash domain, if it is defined by a 
collection of open covers {Ui}fL 0 and associated maps {9±, 9\, 02, • 0l} satisfying the properties 

(l)-(5) above. Because each of the maps is an H 3 dijfeomorphism, we say that the splash domain 
tt s defines a self-intersecting generalized H 3 -domain. 

9.2. An approximating sequence of non self-intersecting domains converging to the 
splash domain. Following “16j, we can then define standard (non self-intersecting) domains fF 
(for e > 0 small enough) by just modifying 9± , and leaving the other charts unchanged. As shown in 
Figure lfT2l our non self-intersecting domain fF will be defined by associated maps {9e±,9 1 ,02, 0l} 

such that 

\\9± - 9±\\ H 3(b+) < Ce, (9.1) 

and such that 

0 < d{6 e + {B + ),Q e _(B + )) < e. (9.2) 


9 + 



FlGURE 9.2. The black dot denotes the point xo where the boundary self-intersects 
(middle). For e > 0, the approximate domain fF does not intersect itself (right). 

In summary, we have approximated the self-intersecting splash domain fl s with a sequence of 
H 3 -class domains fF converging toward 12, such that for each e > 0, <91F does not self-intersect. As 
such, each one of these domains IF, e > 0, will thus be amenable to our local-in-time well-posedness 
theory for free-boundary incompressible Navier-Stokes equations. 
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10. Existence of a splash in finite time in a domain arbitrarily close to a given 

SPLASH DOMAIN 

We next define an initial velocity field of the same type as in Section 16.11 Due to (19.11) , the 
estimates of Section [7] remain unchanged. Similarly, the main proof of Section [ 8 ] works in a similar 
manner due to (19.21) . leading to the necessity of self-intersection at a time T e G (0, lOe). Note 
that since the tangent plane at the intended splash singularity xq is the horizontal plane {. T3 = 0 }, 
d[0-{B + )\ is very close to {*3 = 0} in a small ball B(x 0 , y/e) for e taken sufficiently small; thus, we 
are using the fact that the almost flat portion of 9_(B + ) is very close to {13 = 0} and contains a 
region of diameter at least y/e. 

Furthermore, 

I \v e (o e ±,T e ) - Ms < II v e (o e ±,T e ) - ey 3 + II 0% - Ms 

rpe 

<11 f v e (0±,t) dt\\ 3 +Ce, (10.1) 

Jo 

where we used the estimate (ED in the above inequality (110.11) : hence, from our estimates in Section 

m . 

|| V e {0±,T e ) - 9± || 3 < CV{M 0 )VF + Ce< CV(M 0 )^~e . (10.2) 

This, therefore, shows that the splash-free surface rj e (Cl e , T e ) is at a distance less than CV(Mo)y/e 
from in H 3 . We have then established the following: 

Theorem 10.1. For any given splash domain fl s of class H 3 , there exists a splash domain Q s 
arbitrarily close in H 3 to £} s , and smooth initial data consisting of a non self-intersecting domain 
of class H 3 and a divergence-free velocity field Uq G H 3 (Q e ) satisfying [Def Uq • N e ] x N e = 0 on d£l e , 
such that the flow map r](x,t) solving the Navier-Stokes equations satisfies ry(9fl e ,T*) = Q s . 

That is, in finite time T* > 0, a splash singularity occurs which is very close to a prescribed self- 
intersecting geometry. 
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